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A comprehensive  analysis  is  presented  of  the  nonequilibrium  ablation  of 
carbon  in  an  atmosphere  whose  elemental  composition  consists  of  Oxygen 
Nitrogen,  Carbon  and  an  inert  species.  The  scope  of  the  investigation 
ranges  from  the  low  surface  temperature  regime,  where  the  ablation 
process  is  dominated  by  kinetically  controlled  oxidation  reactions,  to  the 
high  temperature  sublimation  regime  dominated  by  rate -controlled  — •S 


vaporization.  A detailed  examination  is  made  of  the  commonly  employed 
assumption  concerning  the  gas-phase  and  gas-solid  reaction  rates.  In 
particular,  calculations  are  presented  for  the  limiting  assumptions  of  fast 
(equilibrium)  and  slow  (frozen)  homogeneous  and  heterogeneous  reactions. 

The  vaporization  from  a liquid  carbon  melt  layer  is  analyzed  and  calcula- 
tions are  presented  for  the  surface  mass  loss  and  melt  layer  thickness. 
Comparisons  are  made  with  the  available  experimental  data  and  recom- 
mendations made  for  future  experimental  investigation. 
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I.  INTRODUCTION 


The  ablation  of  carbon  may  be  classified  as  equilibrium  or  nonequilib- 
rium, depending  upon  the  nature  of  the  heterogeneous  surface  reactions  and 
the  homogeneous  gas-phase  reactions.  If  the  gas-phase  reactions,  as  well  as 
the  heterogeneous  carbon  vaporization  and  oxidation  reactions,  are  in  equi- 
librium, the  ablation  process  is  termed  equilibrium,  and  the  mass  transfer 
coefficient  is  obtained  from  the  diffusion- controlled  approximation  (e.g.. 

Ref.  1).  The  equilibrium  assumption  is  easily  utilized  and  sufficiently  accu- 
rate for  many  applications.  There  are,  however,  situations  wherein  the 
carbon  ablation  process  is  dominated  by  nonequilibrium  effects,  such  as 
kinetically  controlled  oxidation  and  sublimation  reactions.  Also  included 
among  the  nonequilibrium  effects  is  the  possibility  of  frozen  gas  phase  reac- 
tions and/or  frozen  heterogeneous  reactions  with  the  solid  carbon.  This 
report  is  concerned  with  the  realm  of  the  nonequilibrium  ablation  of  carbon. 

The  available  work  on  nonequilibrium  carbon  ablation  suffers  from 
several  major  and  minor  limitations.  First,  the  scope  of  the  previous 
investigations  is  generally  limited  to  the  oxidation  or  sublimation  regime 
instead  of  both  regions  being  treated  simultaneously. 

A second,  and  more  serious,  criticism  concerns  the  postulated  model 
for  the  gas-phase  and  gas-solid  reactions.  Due  to  a complete  lack  of  kinetics 
data  for  most  of  the  reaction  rates  of  interest,  it  is  necessary  to  invoke  either 
the  frozen  or  equilibrium  assumption  for  many  of  the  heterogeneous  and 
homogeneous  reactions.  Various  combinations  of  the  frozen  and  equilibrium 
assumptions  have  been  utilized  to  obtain  the  ablation  rate.  However,  there 
exists  no  systematic  examination  of  the  influence  of  these  assumptions  upon 
the  calculated  ablation  rates.  Apparently  the  only  investigation  concerning 


Crowell,  P.G.,  "Analytic  Solutions  for  Carbon  Sublimation  in  Atmospheres 
of  Elemental  Carbon,  Hydrogen,  Nitrogen  and  Oxygen,  " Report  No. 
TR-0078(3550-15)-2,  The  Aerospace  Corporation,  El  Segundo,  Calif., 
December  1977. 
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the  influence  of  the  postulated  chemistry  model  is  a limited  study  by 

O .... 

Maurer  to  examine  the  effect  of  the  frozen  versus  equilibrium  assumption 
for  the  gas-phase  carbon  sublimation  species.  The  sensitivity  studies 
contained  in  the  literature  deal  with  the  influence  of  the  allowed  species, 
the  thermochemical  data  used  for  each  species,  the  vaporization  coefficients 
for  the  sublimation  reactions,  et  cetera.  While  such  investigations  are 
necessary,  it  is  felt  that  the  basic  choice  of  a model  for  the  reactions  (frozen/ 
equilibrium)  may  also  exert  significant  influence  on  the  ablation  rate. 

A third  consideration  concerns  the  possible  melting  of  carbon  and  its 
impact  upon  the  predicted  ablation  rate.  The  only  work  on  the  influence  of 

3 

carbon  melting  is  an  equilibrium  study  by  Kratsch  . Although  the  exact 
melt  temperature  is  not  known,  several  experimental  investigations  yield 
values  ranging  from  3800°K  to  4300°K.  Surface  temperatures  in  this  range 
are  often  predicted  utilizing  models  which  assume  ablation  from  solid  carbon 
only.  Furthermore,  there  have  been  suggestions  that  the  "mechanical" 
erosion  observed  in  some  ground  tests  is  actually  due  to  molten  carbon  being 
stripped  off  the  ablating  specimen  by  aerodynamic  forces.  If  molten  carbon 
does  form,  the  ablation  rate  will  be  different  than  that  obtained  from  a solid 
carbon  surface  due  either  to  mechanical  removal  of  the  melt  or  a difference 
in  vapor  pressures  in  going  from  the  gas-solid  system  to  the  gas-liquid 
system.  In  any  event,  this  is  an  area  that  warrants  further  investigation. 

The  fourth  critique  of  the  available  work  involves  the  solution  procedure 
used  to  obtain  numerical  results.  In  many  instances,  the  numerical  method 
is  not  described.  In  those  cases  in  which  it  is  described,  the  method  involves 
a cumbersome,  time-consuming  scheme  with  multiple  iteration  loops.  Most 
of  the  difficulty  experienced  by  some  investigators  in  obtaining  numerical 
results  is  due  to  their  choice  of  chemistry  model  for  the  gas-phase  reactions. 

“^Maurer,  R.E.,  "An  Evaluation  of  Some  Uncertainties  in  Carbon  Sublimation 
Kinetics,  " Aerotherm  TM-73-27,  January  1973. 

^Kratsch,  K.  M.,  "Graphite  Fusion  and  Vaporization,  " ASME/AIAA  10th 
Structures,  Structural  Dynamics  and  Materials  Conference,  April  1969. 

Baker,  R.  L.,  Private  communication 
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That  is,  the  impact  of  the  gas-phase  frozen  versus  equilibrium  assumption  on 

the  numerical  solution  technique  is  apparently  not  recognized.  Thus,  while 
4 5 

Dolton  and  Baker  required  a multiple,  nested  iteration  scheme  for  the 
frozen  assumption,  Ziering^  needed  only  a single  iteration  for  the  equilibrium 
case.  The  solution  procedure  is  not  important  if  the  intent  is  to  generate  a 
set  of  tables  of  mass  transfer  parameter  (B')  for  input  into  a thermodynamic 
response  code.  However,  since  the  mass  transfer  parameter  is  a function  of 
three  variables  (pressure,  temperature,  and  heat  flux),  this  approach  would 
require  a very  tedious  three-dimensional  interpolation  procedure.  A more 
desirable  approach  is  the  utilization  of  a surface  chemistry  subroutine  called 
from  the  thermodynamic  response  code.  In  this  situation,  the  computation 
speed  obviously  becomes  critical. 

Whether  the  minor  limitations  alluded  to  earlier  are  indeed  minor  or 
major  depends  on  the  particular  situation.  A case  in  point  involves  the  solu- 
tion of  the  surface  energy  balance  and  the  utilization  of  the  steady- state 
assumption.  If  the  surface  temperature  is  a specified  boundary  condition, 
then  the  energy  balance  is  not  solved  and  the  steady-state  assumption  is  not 
required.  If,  on  the  other  hand,  the  temperature  must  be  obtained  as  part  of 
the  solution,  the  energy  balance  is  required,  and  the  heat  flux  conducted  into 
the  ablating  surface  must  be  modeled  in  some  manner.  In  order  to  uncouple 
the  solution  of  the  surface  energy  balance  from  the  in-depth  conduction 
problem,  the  steady-state  assumption  is  frequently  used  to  approximate  the 
heat  flux  conducted  into  the  interior.  This  assumption  is  usually  justified  on 
the  basis  of  simplicity  and  reasonable  accuracy.  In  many  situations,  however, 

4 

Dolton,  T.A.,  R.  E.  Maurer,  andH.E.  Goldstein,  "Thermodynamic  Per- 
formance of  Carbon  in  Hyper  thermal  Environments,  " AIAA  3rd  Thermo- 
physics  Conference,  Paper  No.  68-754,  June  1968. 

Baker,  R.  L.,  "Graphite  Sublimation  Chemistry  Nonequilibrium  Effects,  " 
AIAA  Journal,  15(10),  October  1977. 

^Ziering,  M.  B.  and  V.  DiCristina,  "Thermomechanical  Erosion  of  Ablative 
Plastic  Composites,  " AIAA  7th  Thermophysics  Conference,  April  1972, 

Paper  No.  72-299. 
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the  requirements  for  the  validity  of  the  steady-state  assumptions  (i.  e., 
large  time  and  small  thermal  penetration  depth)  are  not  satisfied,  and  its 
use  will  result  in  a considerable  overprediction  of  the  recession  rate. 


Another  "minor”  criticism  of  the  nonequilibrium  ablation  studies  con- 
cerns the  inconsistencies  present  in  many  of  the  formulations.  These  in- 
consistencies arise  as  a consequence  of  the  choice  of  chemical  reaction  rate 
model  (i.e.,  frozen  or  equilibrium)  used  for  the  gas-phase  and  heterogeneous 
reactions.  That  is,  the  use  of  a mixed  set  of  rate-controlled  and  diffusion- 
controlled  heterogeneous  surface  reactions  leads  to  arbitrary  restrictions 
being  imposed  on  the  allowable  gas-phase  reactions  in  order  that  the  rate- 
controlled  reactions  not  be  short-circuited.  Similar  difficulties  arise  when 
the  carbon  sublimation  species  are  frozen  in  the  gas  phase.  Whether  the 
inconsistencies  are  numerically  significant  depends  on  the  particular  set  of 
assumptions  and  allowable  species. 

It  is  the  purpose  of  this  report  to  critically  examine  each  of  the  above 
deficiencies  in  the  existing  formulations  for  the  nonequilibrium  ablation  of 
carbon.  In  particular,  detailed  computations  will  be  carried  out  to  examine 
various  assumptions  concerning  the  chemical  model,  and  comparisons  will  be 
made  with  the  existing  ablation  data.  A computationally  fast  solution  scheme, 
suitable  for  subroutine  use,  will  be  developed  for  the  recommended  formula- 
tion. Recommendations  for  future  development  and  testing  will  be  made  to 
address  those  uncertainties  that  impact  the  predicted  ablation  of  carbon. 

Before  proceeding  to  the  formulation  of  the  next  section,  some  ground 
rules  should  be  established.  The  interest  of  this  report  is  the  ablation  of 
carbon  in  an  atmosphere  whose  elemental  composition  consists  of  carbon, 
nitrogen,  oxygen,  and  an  inert  element  (denoted  by  the  symbol  I).  This 
system  clearly  includes  air  (oxygen,  nitrogen,  and  argon)  as  a special  case. 

If  the  gas  phase  is  assumed  to  be  in  equilibrium,  the  boundary  layer  edge 
elemental  mass  fractions  will  be  presumed  to  be  known.  Furthermore,  the 
edge  species  mass  fraction  will  be  assumed  known  for  any  frozen  species. 
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The  surface  pressure,  temperature,  and  convective  heat  transfer  (either 

blowing  or  nonblowing)  are  considered  known.  Only  the  dominant  gas -phase 

species  (as  determined  from  previous  analyses)  for  a system  of  this  elemental 

composition  will  be  considered  as  well  as  the  first  five  carbon  sublimation 

7 

species  (C^ — *"C^).  Unless  otherwise  indicated,  JANNAF  values  will  be  used 
for  all  species  thermochemical  data.  Since  it  is  intended  that  this  formula- 
tion be  used  in  conjunction  with  a simplified  boundary  layer  analysis,  the  usual 
simplifying  assumptions  (unity  Lewis  and  Prandtl  numbers  and  equal  diffusion 
coefficients)  will  be  invoked. 


! i 

I 


I 


r 


I 


! 


t 


j 


JANNAF  Thermochemical  Tables,  National  Bureau  of  Standards,  NBS-37, 
June  1971. 
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II.  ANALYSIS 


i 


I 


As  indicated  above,  the  major  objective  of  this  report  is  to  examine  the 
influence  of  the  assumed  chemical  state  upon  the  predicted  mass  transfer. 

This  examination  is  carried  out  in  Subsections  II-A  through  II-D.  The  first 
two  subsections  assume  gas -phase  equilibrium  with  frozen  (Subsection  II-A) 
and  equilibrium  (Subsection  II- B)  heterogeneous  carbon-nitrogen  reactions. 

In  Subsection  II-C,  the  carbon  sublimation  species  are  frozen  in  the  gas  phase. 
The  analyses  of  Subsections  II-A  through  II-C  ignore  the  possibility  of  the 
formation  of  a carbon  melt  layer.  The  carbon  melt  layer  problem  is  con- 
sidered in  Subsection  II-D.  Finally,  the  validity  of  the  steady-state  assump- 
tion is  examined  in  Subsection  II- E.  Much  of  the  analysis  is  independent  of 
the  assumed  chemical  state;  this  analysis  is  considered  next. 

The  problem  of  interest  is  to  calculate  the  mass  transfer  from  an 
ablating  carbon  surface  (solid  or  liquid)  over  which  a boundary  layer  is 
flowing.  The  elemental  mass  fractions  at  the  boundary  layer  edge  are  pre- 
sumed known  and  must  sum  to  unity. 


K 


N 


+ K, 


+ K, 


+ K.  = 1 


(1) 


where  subscript  I denotes  the  inert  species. 

At  the  wall  (solid-gas  or  liquid-gas  interface),  the  elemental  mass 
fractions  are  specified  in  terms  of  the  species  of  interest  (see  Appendix  B). 


28  ^CO  44 

w 2 2 

26  52  76  38 

w 2 2 2 4 2 2 

K -i^K  +i^K  +Mk  + 

28  ^CO  44  CO-  26  CN  52  C,N- 
w 2 2 2 


48  . 24 

'4*'2  '^2‘ 


^ 38  ^C.,N 


C. 


Some  comments  about  the  chosen  species  are  in  order.  The  tempera- 
ture range  of  interest  extends  from  540  to  9000°R  and  the  pressure  range 
from  0.  10 — ►BOO  atm.  All  these  species  are  numerically  significant  in 
various  portions  of  the  indicated  temperature  and  pressure  ranges.  Since  a 
unified  treatment  is  sought  over  the  entire  range,  all  these  species  must  be 
considered  simultaneously.  As  will  be  demonstrated  later,  this  set  of  species 
does  not  pose  any  particular  computational  difficulties.  Although  atomic 
oxygen  and  nitrogen  are  not  dominant  species,  they  may  influence  the  abla- 
tion rate  through  kinetically  controlled  surface  reactions  and  should  be 
carried  for  generality.  An  inert  species  is  included  for  generality,  since  it 
is  present  in  many  atmospheres  of  interest.  The  inert  species  does  not 
react  chemically,  since  the  surface  temperature  is  not  high  enough  to  cause 
significant  ionization,  but  its  presence  does  influence  the  molecular  weight. 
Although  the  formulation  contains  only  a single  inert  species,  the  analysis  is 
easily  generalized  for  any  number  of  inert  species.  This  may  be  done  by 
simply  redefining  the  ratio  Kj  /Mj  as  the  summation  of  this  quantity  over  all 
the  inert  species  and  equating  Kj  to  the  summation  of  all  the  inert  elemental 
niass  fractions. 


The  wall  elemental  mass  fractions  may  be  related  to  the  boundary  layer 
edge  values  by  invoking  the  equivalence  between  mass  and  heat  transfer 
(assuming  equal  diffusion  coefficients  and  unity  Prandtl  and  Lewis  numbers). 
This  yields 


I f B' 
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1 + B' 
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P u C 
e e 


H 


(3) 
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The  partial  pressure  of  any  constituent  is  given  as 


I 


i 


I 


PMK. 

P.  = 

i M. 


(4) 


Equation  (2)  is  cast  in  terms  of  partial  pressures  by  multiplying 
through  by  the  product  of  pressure  and  molecular  weight.  This  operation,  in 
conjunction  with  the  substitution  of  Eq.  (3),  yields 

Kq  = f’o  * Pq,  “>  Pco  * Pco, 

e 2 2 

Kn  = Pn  * "8  Pn,  + PcN  * "8  Pc,N,  * 

e 2 2 2 

+ 28  P_  „ +14  P_ 

C4N2  C2N 


PM 


1 + 


:) 


T I B'  + 1 = 12  P^Q  + 12  P^^  + 12  P^„  + 24  P, 


CO. 


CN 


"=2^2 


1 1 


PM  ~ _ 

1 + B'  I 

e 


(5) 


The  sum  of  the  partial  pressures  must  equal  the  mixture  pressure. 


P = P.,  + P.,  + P_.,  + P„  + P_  + P„  K,  + P^  + P^  + 
N N2  CN  C2N2  ^2^  ® ®2 


^co  ^ ^co,  ^ ^^C. 

Z 1 


(6) 


u 
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The  kinetically  controlled  vaporization  reactions  are  modeled  with  the 
Knudsen- Langmuir  equation.®  This  yields  the  mass  loss  due  to  sublimation 


The  heterogeneous  oxidation  reactions,  involving  both  atomic  and 
molecular  oxygen,  are  also  kinetically  controlled  and  are  represented  by 
the  following: 


) 

”°2  ' ""“z  ’’“z) 

In  many  instances,  the  oxidation  reactions  are  represented  as  follows: 

fo  = ('’0°  - J 

/ "o 


where  the  coefficients  K„  and  K_  are  Arhenius  functions,  and  the  exponents 

O Op 

and  n^  indicate  the  order  of  the  reaction. 

O O^ 

Note  that  the  equilibrium  values  of  the  partial  pressures  are  usually 
deleted  from  Eq.  (9).  This  is  because  they  are  usually  several  orders  of 
magnitude  smaller  than  the  actual  partial  pressures.  They  are,  however, 
included  here  to  provide  a formal  means  of  zeroing  the  oxidation  mass  loss 
as  heterogeneous  equilibrium  is  approached. 

^Kennard,  E.H.,  "Kinetic  Theory  of  Gases,  ” McGraw-Hill  Book  Co.,  Inc., 
New  York,  1938. 
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The  equilibrium  values  for  the  partial  pressures  may  be  obtained  from 


f 

k 


Kp  (2C*  + CO)  = 


CO 


eq 


Kp  (C*  + OiZ*CO)  = p-^ 

O 

eq 


(10) 


where  the  pressure  equilibrium  constants  (Kp)  are  defined  for  the  indicated 
reactions. 

In  principle,  the  heterogeneous  carbon  nitridization  reactions  involving 
both  atomic  and  molecular  nitrogen  could  also  be  treated  as  kinetically 
controlled.  That  is,  the  nitridization  mass  loss  could  be  modeled  in  the  same 
manner  as  the  oxidation  through  the  use  of  expressions  analogous  to  those  of 
Eqs.  (8)  through  (10).  Unfortunately,  the  almost  complete  lack  of  kinetics 
data  on  heterogeneous  carbon-nitrogen  reactions  precludes  this  approach.  In 
the  absence  of  reaction  rate  data  for  carbon  nitridization,  recourse  must  be 
made  to  the  limiting  assumptions  of  infinitely  fast  surface  reactions  (i.  e., 
heterogeneous  equilibrium)  or  infinitely  slow  surface  reactions  (i.  e.,  frozen 
heterogeneous  reactions).  It  is  this  choice  of  fast  or  slow  surface  nitridiza- 
tion reactions  which  characterizes  the  earlier  nonequilibrium  ablation  studies. 

4 2 

For  example,  Dolton  and  Maurer  assumes  fast  reactions,  whereas 
Ziering^  assumes  slow  reactions.  As  indicated  above,  the  slow  nitridization 
assumption  is  examined  in  Subsection  II-A,  and  the  fast  reaction  assvimption 
in  Subsection  II- B. 

A literature  search  yielded  only  one  reference  (9)  on  the  experimental 
determination  of  heterogeneous  carbon-nitrogen  reaction  rates.  This  data 
has  been  analyzed  to  determine  which  of  the  two  limiting  nitridization 

^^Goldstein,  H.  W.,  "The  Reaction  of  Active  Nitrogen  with  Graphite,  " 

Journal  of  Physical  Chemistry.  ^(1),  January  1964. 
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aasumptiona  is  most  realistic.  The  test  consisted  of  heating  a graphite 
specimen  to  2365°K  and  impinging  a stream  of  molecuUr  nitrogen  onto  its 
surface  at  a pressure  of  0.  00256  atm.  The  test  was  run  for  30  minutes, 
and  the  gas  phase  was  condensed  out  in  liquid -nitrogen  traps  and  examined 
for  the  presence  of  cyanogen.  A mass  transfer  prediction,  assuming  fast 
heterogeneous  reactions,  indicated  that  the  amount  of  cyanogen  produced  was 
several  orders  of  magnitude  larger  than  the  detection  threshold  of  the 
experimental  setup.  The  calculation  was  repeated,  assuming  slow  hetero- 
geneous reactions,  and  yielded  an  amount  of  cyanogen  about  equal  to  the 
detection  limit.  Since  no  cyanogen  was  detected,  it  is  concluded  that  the 
nitridization  reactions  are  closer  to  the  frozen  limit  than  to  the  equilibrium 
limit.  Due  to  uncertainties  about  the  experiment  and  the  assumptions 
required  to  analyze  the  data  (particularly  the  estimate  for  the  mass  transfer 
coefficient),  this  conclusion  can  hardly  be  viewed  as  definitive.  However, 
in  the  absence  of  any  other  information,  this  result  indicates  that  it  is  more 
reasonable  to  ignore  heterogeneous  nitridization  reactions  than  to  assume 
that  they  are  in  equilibrium. 


Since  the  nitridization  reactions  are  not  assumed  to  be  kinetically 
controlled,  the  total  mass  loss  is  the  sum  of  the  sublimation  and  oxidation 
components.  The  mass  loss  is  expressed  in  terms  of  B'  by  dividing  through 


The  heat  transfer  (p  u C„)  is  obviously  dependent  upon  the  fluid  dynamic 
6 6 n 

aspects  of  the  particular  application  of  interest.  Since  it  is  the  intent  of  this 
report  to  uncouple  the  surface  mass  transfer  problem  from  the  solution  of 
the  boundary  layer  equations,  the  heat  transfer  must  be  modeled  in  some 
manner.  One  approach  is  simply  to  solve  the  mass  tran.afer  problem  for  spe- 
cified values  of  However,  since  it  is  itself  a function  of  the  blowing 

rate,  it  is  more  useful  to  obtain  solutions  in  terms  of  the  nonblowing  heat 
transfer. 
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The  blowing  and  nonblowing  heat  transfer  will  be  related  by 


1 + 


0.60  laminar 


0.  30  turbulent 


In  a number  of  instances,  experimental  data  has  been  interpreted  using 
the  following  result  for  laminar  flow  (ref.  ZZ) 

_ /n(l  -I-  1.28  B')  .... 

C„  ~ 1. 28  B' 


For  values  of  B'  up  to  about  three,  both  these  expressions  yield  virtually 
identical  results  for  laminar  flow.  When  analyzing  data  which  was  reduced 
with  Eq.  (13),  this  report  will  use  the  same  expression.  Otherwise,  the 
expression  of  Eq.  (12)  will  be  used. 

This  completes  the  analysis  common  to  all  the  subsections  that  follow. 
The  specific  reactions  are  peculiar  to  each  subsection  and  are  considered  as 
required. 

A.  EQUIUBRIUM  GAS  PHASE--SIjOW  NITRIDIZATION 


REACTIONS 

In  this  section,  the  mass  transfer  formulation  is  closed  by  assuming 
that  all  gas-phase  reactions  are  in  equilibrium  and  that  heterogeneous  nitridi- 
zatien  reactions  take  place  so  slowly  that  they  may  be  ignored.  Thus, 
carbon-nitrogen  species  are  formed  solely  through  gas-phase  reactions, 
Ziering^  was  apparently  the  first  investigator  to  model  nonequilibrium  car- 
bon ablation  in  this  manner.  This  formulation  differs  from*  Ziering's  in 
that  a larger  number  of  species  is  considered,  and  the  oxidation  regime  is 
treated  simultaneously  with  the  sublimation  regime. 

^^Dorrar.ce,  W.H.,  "Viscous  Hypersonic  Flow,  ■ McGraw-Hill  Book  Co.,  Inc., 
New  York,  1962. 
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The  gas-phase  reactions  are  given  by 


N2  + 2C^2  CN 

CN 

C^N2i=i2  CN  + 2C 

C2NirtCN  + C 

C025=tC  + 20 

C05:±C  + 0 

O2i=!r20 

N2i=?2N 

00 

(14) 

where  the  pressure  equilibrium  constants  are 

for  the  indicated 

reactipiiS. 

The  partial  pressures  for  the  indicated  reactions  are  related  by 


The  pressure  equilibrium  constants  relate  the  partial  pressures  as 
follows: 


P P 
9 C. 


"^11  " ^c' 


P ” P 
10  ^C- 


p p 
^12 


Substitution  of  Eq.  (15)  into  the  oxygen  and  nitrogen  elemental  conser- 
vation equations  of  Eq.  (5)  yields 


PM  7 

= 16  (1  +a^)  Pq  + 32  (05  +0..^)  P^ 


^ J Ivl  ? 

= 28  (a^  +02  +0I3)  P^J^  + 14  (1  +0^  +ag) 


Similarly,  the  carbon  mass  balance  of  Eq.  (5)  becomes 


PM 
1 + B 


,(b. 


12a^Po  + IZa^o  + (12  + 24<>4)Pcn  + 


t (24«^  t 480,)?^^,^  + E“cPc. 


Equation  (6)  yields 


'*  +“6>Po  + '“5  *“7''’o  + <"l  +“2  +“3'PCN  + <1  +“4  + VCN  + 


tV  mt  = p - i:p< 


i 


This  completes  the  problem  formulation.  Specifically,  the  set  of 
equations  to  be  solved  consists  of  Eqs.  (11),  (12)  [or  (13)],  and  Eqs.  (18) 
through  (21).  Since  one  of  the  objectives  is  to  examine  computational  effi- 
ciency, the  solution  procedure  is  now  considered  in  detail. 

A brief  examination  of  the  system  of  equations  makes  it  obvious  that  an 
iterative  solution  procedure  must  be  employed.  Since  the  pressure  equilibrium 
constants  are  known  functions  of  temperature,  it  is  apparent  [see  Eq.  (18)] 
that  guessing  a partial  pressure  for  any  one  of  the  carbon  species  immediately 
yields  the  partial  pressure  values  for  the  remaining  carbon  species.  Since 
Cj  is  usually  the  dominant  carbon  species,  we  adopt  Ziering's  approach  and 
iterate  on  the  partial  pressure  of  C,.  Thus,  for  a guessed  value  of  P , 
the  remaining  are  evaluated  from  Eq.  (18)  and  the  a coefficients 
from  Eq.  (16).  Next,  Eqs.  (19)  and  (21)  are  solved  for  P^,  and 

PM/(1  +B').  Then  Eq.  (20)  is  solved  for  B',  the  blowing  heat  transfer  is 
evaluated  from  Eq.  (12)  [or  (13)],  and  Eq.  (11)  is  also  solved  for  B'.  If  the 
values  of  B'  obtained  from  Eqs.  (11)  and  (20)  agree  to  within  some  prescribed 
tolerance,  the  iteration  has  converged.  If  not,  a new  guess  is  made  for  the 
partial  pressure,  and  the  process  repeated  until  convergence  is  achieved. 
The  Newton -Raph son  iteration  is  used,  and  the  value  of  P-  is  constrained  to 

t>3 

be  less  than  the  upper-bound  value  which  is  the  equilibrium  vapor  pressure  of 
or  the  mixture  pressure,  whichever  is  smaller. 

The  only  computational  problem  encountered  during  the  iteration  is 
the  solution  of  Eqs.  (19)  and  (21).  That  is,  these  three  equations  may  be 
combined  to  yield  a quartic  for  Pq  or  Although  an  analytic  solution 

exists  for  a quartic,  it  was  determined  to  be  computationally  faster  to 
solve  it  iteratively.  This  iteration  proceeds  as  follows:  guess  a value  for 
P Q and  solve  Eq.  (19)  for  PM/(1  + B')  and  P^- The  guess  is  checked  against 
Eq.  (21).  If  the  left-hand  side  is  greater  than  the  right-hand  side  (which  is 
known),  the  guessed  value  of  Pq  is  too  large.  If  it  is  smaller,  the  value  of 
Pq  is  too  small.  This  procedure  continuously  defines  upper-  and  lower- 
bound  values  for  Pq  which  are  used  to  constrain  the  next  Newton- Raphson 
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guess  for  P^.  An  initial  upper-bound  value  for  is  easily  obtained  by 
zeroing  the  CN  partial  pressure  in  Eq.  (21)  and  combining  it  with  the  oxygen 
mass  balance  of  Eq.  (19).  The  solution  is  iterated  until  Eq.  (21)  is  satisfied. 
Generally,  this  iteration  converges  within  three  or  four  guesses. 

A few  closing  comments  are  in  order.  The  only  significant  computa- 
tional difference  between  this  work  and  Ziering's  is  the  necessity  of  solving  a 
quartic,  instead  of  a quadratic,  inside  the  iteration  loop  for  P^  . This  is 
due  to  the  inclusion  of  and  CO^  in  the  allowed  species.  If  on^y  the  high- 
temperature  sublimation  regime  were  of  interest,  these  two  would  be  trace 
species  and  thereby  deleted  from  the  formulation  by  zeroing  and  Ot^.  In 
this  manner,  Eqs.  (19)  and  (21)  would  yield  a quadratic  in  instead  of  a 

quartic.  Similarly,  the  computational  effort  could  be  reduced  by  splitting 
the  sublimation  and  oxidation  regimes  into  separate  high-  and  low-temperature 
domains,  solving  each  domain  separately,  and  then  patching  them  together. 
However,  one  of  the  objectives  of  this  report  is  to  develop  a single  unified 
analysis  valid  for  any  temperature  regime.  It  is  felt  that  the  increased 
generality  is  well  worth  the  slight  computational  penalty  associated  with 
solving  the  quartic. 

Section  I indicated  that  various  approaches  employing  mixed  frozen  and 
equilibrium  assumptions  result  in  inconsistencies  in  the  formulation.  Note 
that  no  such  inconsistencies  are  present  in  the  formulation  of  this  subsection. 
That  is,  the  frozen  nitridization  assumption  does  not  require  the  imposition 
of  any  artificial  constraints  upon  the  equilibrium  gas-phase  reactions. 

This  is  not  the  situation  with  respect  to  the  equilibrium  nitridization  assump- 
tion, as  will  be  seen  in  Subsection  II-B. 

I 
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B.  EQUILIBRIUM  GAS  PHASE- -FAST  NIT RIDIZATION 
REACTIONS 

In  this  subsection,  the  alternate  limiting  assumption  for  the  hetero- 
geneous nitridization  reactions  is  examined.  In  Subsection  II-A,  these 
reactions  were  assumed  to  take  place  so  slowly  that  they  could  be  ignored. 
Here,  it  is  assumed  that  the  reaction  rates  are  sufficiently  fast  to  insure 
heterogeneous  equilibrium.  Specifically,  it  is  assumed  that  the  following 


reactions  are  in 

equilibrium: 

=!= 

2C 

+ N2i=?:2  CN 

* 

❖ 

2C 

+ N^il^C^N^ 

^P 

=!= 

4C 

+ N2^-C4N2 

^P 

>1= 

4C 

+ N^5=i2  C^N 

* 

2 2 

P 

(22) 


where  C denotes  solid  carbon,  and  the  superscript  asterisk  on  the  equilibrium 
constants  indicates  a heterogeneous  reaction.  The  partial  pressures  are 
related  through  the  equilibrium  constants  by 


K. 


K. 


CN 

D 

N, 


^4^2 


K. 


K, 


C^N 


N. 


(23) 
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By  rearrangment  of  these  expressions,  the  f lowing  relationships  are 


obtained: 


p = Of; 

^N2  1 CN 


* 2 

P = Of  P 

C2N2  2 CN 


P “ Of* 

C4N2  " 3 ^CN 


P = Of  P 
C2N  4 CN 


where  the  Of  coefficients  are  defined  by 


“1  " 

^P, 


* 

«4  " 


With  the  exception  of  the  homogeneous  carbon-nitrogen  reactions,  the 
gas-phase  reactions  are  identical  to  those  of  Subsection  II-A.  However,  the 
first  four  reactions  of  Eq.  (14)  must  be  deleted  in  order  to  avoid  short- 
circuiting  the  nonequilibrium  sublimation  reactions.  That  is,  inclusion  of 
these  homogeneous  carbon-nitrogen  reactions,  together  with  the  heterogeneous 
ri-actions  of  Eq.  (22),  would  force  the  partial  pressure  of  each  carbon  subli- 
mation sjjecies  to  be  equal  to  its  corresponding  vapor  pressure,  thereby 
zeroing  the  Knudsen- Langmuir  equations  and  imposing  equilibrium  upon  the 
sublimation  reactions.  In  order  to  prevent  this,  all  gas-phase  reactions 


1 
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involving  the  carbon  sublimation  species  and  nitrogen  or  carbon-nitrogen 
species  must  be  suppressed  (i.e.,  assumed  to  be  frozen).  For  example, 
the  following  reactions  are  not  allowed  to  be  in  equilibrium: 

C + N 

C + CN5=tC2N 

(26) 

The  only  gas -phase  reactions  involving  the  carbon  sublimation  species  which 
are  allowed  to  be  in  equilibrium  are  the  dissociation  reactions  of  Eq.  (17) 
and  reactions  involving  carbon-oxygen  compounds. 

These  restrictions  on  the  allowable  gas -phase  reactions  are  accomplished 
simply  by  replacing  the  first  four  homogeneous  reactions  of  Eq.  (14)  by  the 
heterogeneous  reactions  of  Eq.  (22).  This  results  in  the  coefficients  Ot, 
through  Ct^  being  replaced  by  Otj  through  01^.  This  is  the  only  change  required 
in  the  analysis  of  Subsection  II-A  to  implement  the  assumption  of  fast  nitridi- 
zation  reactions. 

It  should  be  pointed  out  that  the  use  of  Eqs.  (22)  does  not  imply  that  the 
carbon-nitrogen  species  are  frozen  in  the  gas  phase.  While  they  are  pro- 
hibited from  reacting  with  the  carbon  sublimation  species,  they  are  allowed 
to  be  in  equilibrium  with  nitrogen  and  with  one  another.  That  is,  Eqs.  (22) 
are  exactly  equivalent  to  any  one  of  the  heterogeneous  reactions  plus  three 
equilibrium  gas -phase  reactions.  For  example,  Eqs,  (22)  could  be  replaced 
by  the  following  set  of  equivalent  reactions: 

2C*  + 2 CN 

C2N25=?^2  CN 

C4N25=!:2  C^N 

^2  <27) 
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In  summary,  the  use  of  the  fast  nitridization  assumption,  in  conjunction 
with  nonequilibrixim  sublimation  reactions,  precludes  the  use  of  a self- 
consistent  set  of  equilibrium  reactions  for  the  gas  phase.  Therefore,  the 
validity  of  the  fast  nitridization  assumption  rests  on  the  existence  of  an 
arbitrary  mixed  set  of  frozen  and  equilibrium  gas-phase  reactions.  This 
seems  somewhat  implausible  and  must  be  regarded  as  casting  additional 
doubt  upon  the  fast  nitridization  assumption. 

C.  FROZEN  CARBON  SPECIES 

In  this  subsection  we  examine  the  assumption  that  the  carbon  sublima- 
tion species  are  frozen  in  the  gas  phase.  As  pointed  out  by  Ziering  , this 
assumption  seems  questionable  in  view  of  the  high  pressures  and  tempera- 
tures of  interest.  The  frozen  assumption  appears  to  stem  more  from 

4 

historical  precedent  (i.e.,  Dolton  ) than  from  any  physical  rationale. 

The  previous  analysis  [Eqs.  (1)  through  (13)]  remains  valid  for  this 
subsection,  but  the  reactions  are  recast  in  somewhat  different  form.  The 
carbon-oxygen  species  and  the  oxygen  and  nitrogen  dissociation  reactions  are 
assumed  to  be  in  equilibrium 


CO  + 


O 


2 


2 O 


N 


2 


2 N 


(28) 


Note  that  the  carbon  monoxide  dissociation  reaction  has  been  deleted  in 
keeping  with  the  requirement  that  gas -phase  reactions  involving  free  carbon 
are  prohibited  by  the  frozen  assumption. 
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The  equilibrium  relationships  are  given  by 


P = <*  P P 
CO^  13^0^00 


P = Ot  p 
^>2  7^0 


P = 0(  P 

8^CN 


where  the  01  coefficients  are  defined  by 

_ - 1 


(29) 


«13  - Kp 


13 


7 K, 


-8 


= Kp 


(30) 


The  carbon-nitrogen  species  are  assumed  to  be  in  gas -phase  equilibrium 
with  each  other.  This  yields 


CN 

C 2 C-N 

H c L 

Since  there  are  four  carbon-nitrogen  species  present  in  the  system,  a 
fourth  reaction  is  required  to  close  the  formulation.  There  is  no  way  to  write 
a fourth  linearly  independent  expression  involving  the  species  of  Eq.  (31). 
Therefore,  the  required  expression  must  involve  any  of  those  species  and 
carbon.  Since  gas -phase  reactions  involving  the  carbon  sublimation  species 
are  prohibited,  the  fourth  equation  must  be  a heterogeneous  nitridization 
reaction.  Therefore,  we  require  the  following  reaction  to  be  in  equilibrium: 

2C’^  + CN  (32) 
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As  indica.ted  in  Subsection  II-B,  Ec^s,  (31)  and  (3Z)  are  exactly  equiva- 
lent to  the  full  set  of  heterogeneous  reactions  of  Eqs.  (22).  Therefore,  we 
adopt  Eqs.  (23)  through  (25)  for  the  formulation  of  this  subsection. 

The  carbon  dissociation  reactions  of  Eq.  (17)  are  deleted  and  replaced 
by  the  following  expressions: 


me 

^e  e^H  ‘-'i  ^i 

w e 


i=l-5 


(33) 


The  frozen  species  equations  are  obtained  from  the  species  continuity 
equation  and  the  analogy  between  mass  and  heat  transfer  (for  equal  diffusion 
coefficients  and  unity  Lewis  and  Prandtl  numbers). 

Casting  Eq.  (33)  in  terms  of  partial  pressures  and  combining  with  the 
Knudsen- Langmuir  Eq.  (7)  yields 


PM 
1 + B' 


P u C„  M„  P„ 
e e H C.  C. 

1 1 


(34) 


M 


C. 


C.  V ZttRT  \^^C.  ■ ^C.J  ^ ^e'^e^H  ^C. 


Combining  the  oxygen  mass  balance  of  Eq.  (5)  and  the  summation  of 
partial  pressures  [Eq.  (6)]  and  using  Eqs.  (24),  (29),  and  (30)  yields 


APq  + BPq  + C = 0 


^C  =2X.P 


^ PM  \ °e 

TTW]  TT 


A = 16  - Oj^)  > 0 


B 


“l3  KqJ  ^ » 


(35) 
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>>  • p - EPc.  - (1  + + «8)PcN  - '“I  + <*I  + PcN  * 


PM  € 
1 + B‘  M, 


(35  cont.  ) 


The  nitrogen  mass  balance  of  Eq.  (5)  becomes 


28(at  +a*  +0*)  + 14(1  + a*  + ajp 


The  carbon  mass  balance  of  Eq.  (5)  is 


8'"CN  1 


PM  - 
L + B* 


(b'  + KcJ  = ‘2(‘  * “i3Po)Pco  + + 2''“>cn  * 

+ (240^’'  + 48a3)P^j^  + Pj.  (37) 

An  examination  of  this  system  of  equations  indicates  that  the  carbon 
species  partial  pressures  cannot  be  evaluated  solely  from  a guessed  value  of 
one  of  them.  This  is  a consequence  of  replacing  the  equilibrium  expressions 
of  Eq.  (17)  by  the  frozen  expressions  of  Eq.  (33).  This  considerably  compli- 
cates the  numerical  solution  procedure,  requiring  a double  iteration  involving 
the  blowing  heat  flux  and  one  of  the  carbon  partial  pressures.  The  outer 

iteration  loop  consists  of  guessing  a value  for  p u C,„  solving  the  inner 

e e H 

iteration  loop  for  B'  and  checking  the  guess  against  Eq.  (12). 

For  a supplied  value  of  p^u^Cj^,  the  iteration  procedure  for  the  inner 
loop  is  described  by  the  following:  Guess  a value  for  P and  solve  Eq.  (34) 

C3 

for  PM/(1  +B'),  Since  Eq.  (34)  is  valid  for  every  sublimation  species,  it 
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may  be  resolved  for  the  remaining  Pp  using  the  value  of  PM/(1  +B')  obtained 

^i 


from  the  guess  for  P 


Next,  Eq.  (36)  is  solved  for  and  Eq.  (35)  for 


Pq  and  Pqq«  The  mass  transfer  parameter  B'  is  then  calculated  from 

Eq.  (37)  and  checked  against  the  value  obtained  from  Eq.  (H).  If  the  two 

values  of  B'  do  not  agree,  a new  guess  is  made  for  the  partial  pressure, 

and  the  process  repeated  until  the  iteration  converges.  The  guessed  value  of 

P^  is  constrained  by  requiring  it  to  be  less  than  the  vapor  pressure, 

an  "d^by  the  requirement  that  the  B and  C coefficients  of  Eq.  (35)  have  the 

indicated  signs.  If  B is  negative,  the  guessed  value  for  P„  is  too  large, 

^3 

whereas  it  is  too  small  if  C is  positive. 

D.  VAPORIZATION  OF  MOLTEN  CARBON 

As  indicated  previously,  the  major  objective  of  this  report  is  to  examine 
the  impact  of  various  assumptions  concerning  the  chemical  state  of  the  sys- 
tem upon  the  predicted  ablation.  Thus,  in  Subsections  II-A  through  II-C,  the 
influence  of  the  limiting  reaction  rates  (i.e.,  fast-slow)  was  examined.  How- 
ever, in  all  cases,  it  was  assumed  that  carbon  sublimed  from  a solid  surface. 
In  this  subsection,  we  shall  examine  the  influence  of  the  surface  state  (i.e., 
solid  or  liquid  phase)  on  the  ablation  formulation.  In  other  words,  it  will  be 
assumed  that  the  carbon  vaporization  reactions  occur  from  a liquid  carbon 
surface. 

3 

Kratsch  used  an  equilibrium  analysis  to  calculate  the  liquid  layer 
mass  loss  and  obtained  values  about  an  order  of  magnitude  larger  than  the 
cc)mpa rable  recession  due  to  sublimation  from  a solid  surface.  It  is  not 
clear  exactly  how  these  large  recession  rates  were  obtained,  but  one  possi- 
bility might  be  that  the  melt  layer  was  mechanically  removed  from  the  sur- 
face by  aerodynamic  forces.  Since  recession  rates  of  this  magnitude  are 
not  indicated  by  experimental  data,  some  doubt  must  be  cast  on  the  validity 
of  the  Kratsch  analysis.  The  lower-bound  limiting  assumption  is  that  the 
melt  layer  remains  in  place,  and  the  mass  loss  is  simply  the  vaporization 
rate  from  the  liquid  surface.  If  the  melt  layer  is  sufficiently  thin,  then 


r 


viscous  forces  and  surface  adhesion  will  resist  mechanical  removal  and 
flowing  of  the  liquid  layer.  This  report  will  assume  that  the  liquid  layer  is 
thin  enough  such  that  it  does  not  flow  and  is  not  removed  by  aerodynamic 
forces. 

It  should  be  apparent,  therefore,  that  the  analysis  of  the  vaporization 
from  a liquid  surface  is  similar  to  the  sublimation  analysis  from  a solid 
surface.  The  only  difference  is  that  the  properties  of  the  solid  carbon  are 
replaced  by  the  properties  of  liquid  carbon.  In  particular,  the  equilibrium 
vapor  pressure  existing  over  a liquid  surface  is  different  than  that  existing 
over  a solid  surface.  Thus,  the  vapor  pressure  for  solid  carbon  must  be 
replaced  by  the  liquid  value.  The  vapor  pressure  for  both  phases  is  obtained 
from  Appendix  A. 

The  phase  state  of  the  surface  probably  influences  other  aspects  of  the 
ablation,  such  as  the  oxidation  and  nitridization  reaction  rates  and  the  values 
for  the  vaporization  coefficients.  Since  the  partial  pressures  of  atomic  and 
molecular  oxygen  are  vanishingly  small  for  the  melt  temperatures  of  interest, 
the  oxidation  reactions  are  insignificant.  The  situation  with  respect  to  the 
vaporization  coefficients  and  nitridization  reactions  is  not  as  obvious.  In  the 
absence  of  any  other  information,  the  same  set  of  vaporization  coefficients 
will  be  used  for  both  the  solid  and  liquid  phases.  The  use  of  the  equilibrium 
assumption  for  the  heterogeneous  nitridization  reactions  yields  the  same 
inconsistencies  as  were  encountered  with  the  solid  carbon  surface.  There- 
fore, the  slow  nitridization  assumption  is  utilized. 

In  summary,  the  analysis  for  the  vaporization  of  a carbon  melt  layer  is 
identical  to  that  of  Subsection  II-A,  with  the  correct  vapor  pressure  law  for 
liquid-gas  equilibrium. 
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E. 


SURFACE  ENERGY  BALANCE 


The  analysis  of  the  preceding  subsections  treated  the  surface  tempera- 
ture as  a known  quantity.  Although  a prescribed  temperature  is  a common 
boundary  condition  for  boundary  layer  calculations,  most  ablation  analyses 
require  the  determination  of  the  surface  temperature  as  part  of  the  solution. 
The  temperature  is  obtained  by  satisfying  the  surface  energy  balance  which 
equates  the  net  aerodynamic  surface  heat  flux  to  the  energy  conducted  into 
the  interior  of  the  body. 

The  surface  energy  balance  may  be  expressed  as 

q„  = p u C„  / 7^  \ Fho  - H - B’  (H  - H„)l  - (TCT'^ 

^C  '^e  e Hq  yCpj  J \_  R w w SJ  w 

H = EKh.  (38) 

where  q^  denotes  the  energy  conducted  into  the  body  and  Hg  is  the  enthalpy  of 
the  surface  phase  (liquid  or  solid)  supplied  by  the  Appendix. 

In  principle,  a knowledge  of  the  interior  temperature  field  is  required  to 
evaluate  q^,  thereby  coupling  the  in-depth  material  response  calculation  to 
the  surface  mass  transfer  calculation.  If  the  conduction  heat  transfer  can  be 
expressed  solely  in  terms  of  the  surface  temperature,  the  mass  transfer  cal- 
culation may  be  uncoupled  from  the  in-depth  material  response,  resulting  in 
a considerable  simplification.  This  uncoupling  is  usually  accomplished  by 
invoking  the  steady-state  assumption  to  evaluate  the  conduction  heat  transfer. 
It  is  the  purpose  of  this  subsection  to  present  an  alternate  model  for  the 
conduction  heat  transfer  which  includes  the  influence  of  the  time  scale. 

It  is  assumed  that  the  conduction  may  be  represented  as  a one-dimen- 
sional process  normal  to  the  surface.  It  is  also  assumed  that  the  body 
dimension  normal  to  the  surface  is  much  larger  than  the  thermal  penetration 
depth  for  the  time  scale  of  interest,  such  that  the  body  may  be  represented 
as  a semi-infinite  domain. 
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The  conduction  equation  may  be  written  in  a coordinate  system  which 
translates  with  the  surface  as  it  ablates.  This  yields 


dT 


dT 


(k£) 


(39) 


Integrating  from  x = 0 (surface)  to  x = oo  and  assuming  constant  material 
properties  yield 


•( 


■ <)... 
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dT 

dt 


dx 


(40) 


where  Tg  is  the  bulk  temperature  of  the  in-depth  material. 

The  first  term  of  this  expression  is  recognized  as  the  steady- state 
approximation.  It  is  apparent  that  the  steady -state  assumption  will  be  valid 
if  the  time  scale  is  large  enough  so  that  the  integral  term  is  negligible. 
However,  if  the  time  scale  is  small  or  the  recession  rate  (S)  is  low,  then  the 
integral  term  dominates,  and  the  steady-state  assumption  is  quite  poor.  In 
order  to  introduce  the  time  scale  into  the  problem,  the  integral  term  must  be 
approximated  in  some  manner.  This  is  accomplished  by  using  the  solution  of 
Carslaw  and  Jaeger  for  a semi-infinite  domain  (of  a single  phase)  with 
constant  surface  temperature  and  recession  rate.  This  yields 
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Carslaw,  H.S.  andJ.C.  Jaeger,  Conduction  of  Heat  in  Solids.  Oxford 
University  Press,  1959,  p.  388. 
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This  expression  is  plotted  in  Fig.  1 where  it  is  apparent  that  the 
requirement  for  the  steady-state  approximation  to  be  valid  is  given  by 

1 v/T  > I (42) 

2 c 

Note  that  Cq.  (41)  possesses  the  correct  limiting  behavior  as  the  time 
and/or  recession  rate  vanishes,  whereas  the  steady-state  model  is  clearly 
incorrect. 

Thus,  for  a specified  value  of  time,  Eq,  (41)  may  be  solved  with  the 
surface  energy  balance  to  yield  the  temperature  and  recession  rate.  The 
time  scale  is  simply  chosen  to  be  the  elapsed  time  from  the  initiation  of  the 

heating  pulse.  It  is  felt  that  the  introduction  of  the  time  scale  in  this  manner  i 

should  substantially  increase  the  accuracy  for  applications  in  which  the 
steady-state  approximation  is  invalid,  while  retaining  the  computational 
simplicity  associated  with  uncoupling  the  in-depth  response  from  the  surface 
mass  transfer. 

Equations  (38)  through  (42)  provide  all  the  necessary  information  to 
analyze  the  ablation  from  a solid  carbon  surface.  If,  however,  the  solid 
carbon  changes  phase  by  melting,  additional  information  is  required.  The 
transient  solution  required  to  obtain  Eq.  (41)  is  only  valid  for  a solid  carbon 
body  and  is  not  applicable  for  the  carbon  melt  layer  problem.  Unfortunately, 
the  analogous  transient  conduction  solution  for  a two-phase  semi-infinite 
body,  with  a receding  surface,  does  not  exist.  Therefore,  the  analysis  of 

the  melt  layer  problem  will  utilize  the  steady-state  conduction  assumption.  ; 

The  expression  of  Eq.  (40)  is  also  valid  for  the  situation  where  a phase 
change  occurs  in  depth.  Thus,  the  steady-state  conduction  assumption  (for 
either  a single  phase  or  two  phases)  may  be  expressed  in  terms  of  enthalpies 

a s I 

1 

I 

(Hg  - Hg)  (43) 

where  is  the  enthalpy  of  solid  carbon  at  the  bulk  temperature.  j 

B . j 
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The  thickness  of  the  liquid  carbon  layer  is  also  required  to  check  the 
assumption  that  the  melt  layer  is  thin.  The  steady-state  conduction  solution 
yields  the  following  result  for  the  melt  layer  thickness  (i.e.,  6) 


in 


(44) 


where  denotes  the  enthalpy  of  liquid  carbon  at  the  indicated  temperature, 
and  is  the  liquid  thermal  diffusivity. 
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III.  RESULTS  AND  DISCUSSION 


The  assumptions  utilized  in  the  analysis  are  examined  in  this  section, 
the  various  models  are  compared  to  each  other,  and  comparisons  are  made 
with  experimental  data.  In  Subsection  lU-A,  calculations  are  presented  for 
the  kinetically  controlled  oxidation  regime,  and  the  impact  of  uncoupling  the 
boundary  layer  calculation  from  the  surface  mass  transfer  calculation  is 
examined.  In  Subsection  III-B,  calculations  utilizing  the  slow  nitridization 
assumption  are  presented.  Comparisons  are  made  with  the  fast  nitridization 
assumption  in  Subsection  III-C  and  the  frozen  sublimation  species  in  III-D. 
Carbon  melt  layer  calculations  are  shown  in  Subsection  III-E,  and  the  influence 
of  the  steady-state  conduction  assumption  is  examined  in  Subsection  III-F. 
Finally,  some  comparisons  with  experimental  ablation  data  are  presented  in 
Subsection  III-G. 

Before  proceeding,  some  ground  rules  for  the  computations  must  be 
established.  First,  all  the  calculations  of  this  report  are  for  air  environ- 
ments. Thus,  the  boundary  layer  edge  mass  fractions  for  the  elements  and 
the  carbon  species  (required  for  the  frozen  assumption)  are  given  by 

Kq  = 0.  232  = 0.  768  Kj  = 0 

® e ^e 

= 0 = 0 (45) 


Although  the  one  percent  mass  fraction  of  argon  present  in  the  atmos- 
phere could  be  included  as  the  inert  species,  its  effect  is  negligible  and  is 
usually  ignored. 
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As  indicated  earlier,  it  is  not  the  intention  of  this  report  to  examine  the 
influence  of  species  property  variations  upon  the  predicted  mass  transfer. 
Therefore,  a fixed  set  of  species  data  was  used  throughout.  With  one  excep- 
tion, the  JANNAF  free  energy  functions  and  nominal  heats  of  formation  were 
utilized.  However,  the  JANNAF  values  for  the  molecule  are  generally 
believed  to  be  in  error.  Accordingly,  the  recent  calculation  of  Hansen  and 
Pearson^^,  as  curve  fit  by  Pearson  and  Davy  , were  used  for  the  free 

energy  function.  Note  that  these  results  are  essentially  identical  to  the 

.14  ...  15 

earlier  values  obtained  by  Straus  and  Thiele  . The  Wachi  and  Gilmartin 
heat  of  formation  was  adopted  for  A total  of  five  carbon  sublimation 

species  was  included  in  the  calculations,  and  the  JANNAF  data  was  used  for 
the  other  four  species. 

According  to  the  JANNAF  data,  the  only  major  uncertainty  in  the  species 
of  interest  is  the  C2^'  heat  of  formation.  That  is,  JANNAF  quotes  a nominal 
heat  of  formation  of  131.  9 kcal/mol,  with  a possible  variation  of  *30  kcal/mol. 
This  represents  a tremendous  spread  in  the  predicted  C2N  mass  fraction; 
it  varies  from  being  the  dominant  species  in  the  system  (for  the  smallest  heat 
or  formation)  to  only  a trace  species  (for  the  largest  heat  of  formation).  A 
more  recent  value  for  the  C,N  heat  of  formation  was  given  by  Meyer  and 
Setser  as  184  kcal/mol.  It  is  clear  that  additional  work  is  required  to 
narrow  this  range  of  values.  In  any  event,  the  nominal  JANNAF  value  was 
used  in  this  report. 

Hansen,  C.  F.  and  W.E.  Pearson,  "A  Quantum  Model  for  Bending  Vibra- 
tions and  Thermodynamic  Properties  of  C3,  " Canadian  Journal  of  Physics, 

1973. 

^^Pearson,  W.E.  and  W.  C.  Davy,  "New  Thermodynamic  Functions  for  the 
C3  Molecule,  " AIAA  Journal,  n(8),  August  1973. 

^*^Strauss,  11.  E.  and  E.  Thiele,  "Thermodynamics  of  C3  General  Methods 
for  Nonrigid  Molecules  at  High  Temperature,  " Journal  of  Chemical  Physics, 
-^(7),  April  1967. 

^^Wachi,  F.  M.  and  D.  E.  Gilmartin,  "Heat  of  Formation  and  Entropy  of  C^ 
Molecule,"  Report  No.  TR -0  1 72  (2250  - 40)  - 7,  The  Aerospace  Corporation, 

El  Segundo,  Calif.,  May  1972. 

Meyer,  J.A.  and  D.  W.  Setser,  "Chemiluminescence  from  the  Reaction  of 
Oxygen  Atoms  with  Dicyanoacetylene,  " Journal  of  Physical  Chemistry,  74, 
1970,  p.  3452. 
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There  remains  considerable  uncertainty  about  the  vaporization  coeffi- 
cients for  the  carbon  species,  with  some  investigators  claiming  values  close 

4 

to  unity  and  others  (see  Dolton  ) obtaining  values  small  compared  to  unity. 

For  unity  vaporization  coefficients,  the  nonequilibrium  formulation  yields 

5 

results  quite  similar  to  the  equilibrium  values  (see  Baker  ).  The  vaporiza- 
tion coefficients  used  by  Dolton  yield  significant  departures  from  equilibrium 
and  are  the  values  used  herein.  These  values  are  given  as 

= 0.24  yS  = 0.50  = 0.023 

^2  '"3 

= 0.25  = 0.0019  (46) 

Note  that  there  is  no  reason  to  believe  that  the  vaporization  coefficients 
are  independent  of  temperature.  They  are  invariably  assumed  constant  due  to 
a lack  of  data  on  the  influence  of  temperature.  Similarly,  the  vaporization 
coefficients  from  a liquid  carbon  surface  would  not  be  expected  to  be  the  same 
as  those  from  a solid  surface.  It  has  been  postulated  that  vaporization  from  a 
liquid  surface  should  take  place  more  readily  than  from  a solid  surface  and, 
therefore,  that  the  liquid  vaporization  coefficients  should  tend  towards  unity. 
However,  Hirth  suggests  that  liquids  with  polymeric  vapor  species  should 
have  values  less  than  unity.  In  the  absence  of  any  other  information,  it  is 
assumed  that  the  Dolton  vaporization  coefficients  also  apply  to  liquid  carbon. 


A.  KINETICALLY  CONTROLLED  OXIDATION 

The  objectives  of  this  subsection  are  twofold.  First,  to  determine  the 
applicability  of  the  iteration  scheme  of  Subsection  II- A at  low  temperatures 
and  second,  to  compare  the  results  with  an  exact  boundary  layer  calculation. 
The  choice  of  the  partial  pressure  as  the  iteration  parameter  was  based 

on  the  fact  that  it  is  the  dominant  carbon  species  at  elevated  temperatures. 

__ 

Hirth,  J.  P.,  "Evaporation  and  Sublimation  Mechanisms,  " The  Characteri- 
zation of  High  Temperature  Vapors,  J.  L.  Margrave  (ed.  ),  John  Wiley  & 
Sons,  Inc.,  New  York,  1967. 
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However,  in  the  low-temperature  oxidation  regime  it  exists  only  as  a trace 
species.  In  some  instances,  the  use  of  a trace  species  as  an  iteration  para- 
meter can  lead  to  very  slow  convergence  or  nonconvergence.  No  difficulties 
of  this  sort  were  encountered  in  the  calculations  reported  herein. 

As  indicated  in  Section  II,  several  assumptions  were  used  to  uncouple 
the  boundary  layer  calculation  from  the  surface  mass  transfer  calculation. 
Specifically,  the  assumption  of  unity  Prandtl  and  Lewis  numbers  and  equal 
diffusion  coefficients  were  employed  to  relate  the  wall  and  edge  elemental 
mass  fractions  [Eq.  (3)].  In  order  to  assess  the  impact  of  these  assumptions, 
comparison  is  made  with  an  exact  boundary  layer  calculation  which  does 
not  utilize  such  assumptions. 

The  case  of  interest  consists  of  a stagnation  point  calculation  at 
moderate  pressure  over  a range  of  temperature  up  to  the  CO  diffusion  limit. 
The  oxidation  reaction  involves  only  molecular  oxygen  and  is  characterized 
by  two  reaction  rates,  designated  as  slow  and  fast.  The  reaction  rates  and 
reaction  order  are  given  by 


II 

-42300/RT 

Slow 

reaction: 

10 

e 

= 6. 73  X 

-44000/RT 

Fast 

reaction: 

10 

e 

Reaction  order: 

"°2 

= 0.  50 

The  values  of  u C,,  required  for  the  formulation  of  Subsection  II-A 
'^e  e H 

were  obtained  froni  the  exact  solution.  Thus,  the  blowing  model  [Eqs,  (12) 
or  (13)]  does  not  enter  into  the  comparison. 


^^Bernstein.  H.  and.T.R.  Baron,  "Stability  Characterization  of  Refractory 
Materials  under  High  Velocity  Atmospheric  Flight  Conditions,  Part  IV, 
Volume  II,  Calculation  of  the  General  Surface  Reaction  Problem,  " AFML 
TR-69-84,  December  1969. 
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The  comparison  is  shown  in  Fig.  2 where  it  is  apparent  that  the  exact 
and  approximate  solution  differ  by  less  than  10  percent.  It  is  concluded, 
therefore,  that  the  assumptions  used  to  uncouple  the  boundary  layer  and  mass 
transfer  calculations  are  valid.  In  fact,  the  uncertainties  in  the  oxidation 
kinetics  will  swamp  any  of  the  uncoupling  assumptions. 

B.  SLOW  NITRIDIZATION  CALCULATIONS 

In  the  previous  subsection,  the  computation  scheme  of  Subsection  Il-A 
was  used  to  obtain  results  in  the  low-temperature  oxidation  regime.  Since 
carbon-nitrogen  compounds  and  carbon  sublimation  species  are  insignificant 
at  low  temperatures,  the  impact  of  the  assumptions  of  Subsections  II-A 
through  II-C  are  unrecognizable  in  the  oxidation  regime.  In  order  to  dis- 
tinguish between  the  various  reaction  rate  assumptions,  a set  of  calculations 
is  presented  for  the  sublimation  regime.  Since  the  atomic  and  molecular 
oxygen  partial  pressures  vanish  at  high  temperatures,  the  choice  of  reaction 
rate  for  the  oxidation  reaction  exerts  negligible  influence  on  the  results. 

The  mass  transfer  parameter  is  presented  in  Figs.  3 through  5 for 
three  pressures  and  a range  of  nonblowing  heat  transfer  coefficients.  The 
calculations  used  the  laminar  blowing  expression  of  Eq.  (12)  and  the  Dolton 
vaporization  coefficients.  For  comparative  purposes,  the  equilibrium  values 
are  also  shown.  As  expected,  as  the  pressure  increases,  larger  values  of 
heat  flux  are  required  to  obtain  significant  departures  from  the  equilibrium 
values . 

C.  COMPARISON  OF  FAST  AND  SLOW  NITRIDIZATION 

ASSUMPTIONS 

The  B'  calculations  for  the  fast  and  slow  nitridization  assumptions  are 
compared  in  Figs.  6 through  9 for  several  pressures  and  nonblowing  convec- 
tive heat  fluxes.  The  laminar  blowing  expression  of  Eq.  (12)  and  the  Dolton 
vaporization  coefficients  were  utilized  for  the  calculations,  and  the  gas-phase 
reactions  were  in  equilibrium,  subject  to  the  restrictions  required  for  the 
fast  nitridization  assumption. 
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The  term  "slow  nitridization,  " in  the  context  used  herein,  means 
infinitely  slow  or,  equivalently,  that  no  mass  loss  occurs  due  to  heterogeneous 
carbon-nitrogen  reactions.  Therefore,  the  fast  nitridization  assumption, 
which  maximizes  the  mass  loss  due  to  heterogeneous  carbon-nitrogen  reac- 
tions, always  yields  a larger  B'  value  than  the  slow  assumption,  for  given 
values  of  temperature,  pressure  and  convective  heat  flux.  It  is  apparent  that 
significant  differences  exist  only  in  a narrow  temperature  range  just  prior  to 
the  onset  of  the  sublimation  reactions.  As  expected,  the  difference  between 
the  fast  and  slow  assumptions  increases  with  increasing  departure  from 
equilibrium  (low  pressures  and  large  convective  heat  fluxes). 

The  maximum  difference  in  B'  for  the  two  nitridization  assumptions  is 
shown  in  Fig.  10  as  a function  of  pressure  and  nonblowirg  heat  flux.  The 
maximum  value  of  the  B'  ratio  is  obtained  from  results  such  as  those  in 
Figs.  6 through  9 at  whatever  temperature  it  occurs.  Substantial  differences 
are  shown  to  occur  for  the  fast  and  slow  assumptions  for  sufficiently  large 
departures  from  equilibrium.  However,  for  the  small  degree  of  nonequilibrium 
experienced  in  ground  test  facilities  and  nominal  ICBM  reentry  flight  tests, 
the  differences  in  mass  transfer  parameter  for  the  two  limiting  assumptions 
are  quite  small.  That  is,  solving  the  surface  energy  balance  for  both  the  fast 
and  slow  nitridization  assumptions  yields  maximum  differences  in  recession 
rate  of  the  order  of  10  percent  for  the  range  of  heat  flux,  pressure  and  recovery 
enthalpy  of  interest.  The  resultant  difference  in  surface  temperature  is  of  the 
order  of  100°R.  Thus,  differences  in  mass  loss  due  to  temperature-dependent 
erosion  mechanisms  will  also  be  negligible. 

D.  FROZEN  CARBON  SUBLIMATION  SPECIES 

Mass  transfer  calculations  are  presented  here  for  the  assumption  that 
the  carbon  species  which  sublime  from  the  surface  are  frozen  in  the  gas  phase. 
As  indicated  in  Subsection  II-C,  this  frozen  assumption  also  requires  utiliza- 
tion  of  the  fast  nitridization  assumption  for  the  heterogeneous  carbon-nitrogen 
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reactions.  Therefore,  the  influence  of  the  frozen  gas -phase  reaction  is  seen 
by  comparison  with  the  equilibrium  gas -phase  calculations  with  the  fast 
nitridization  assumption.  This  comparison  is  presented  in  Figs.  11  through 
13  as  a function  of  temperature,  pressure,  and  nonblowing  heat  transfer 
coefficient. 

The  difference  between  the  frozen  and  equilibrium  assumptions  for  the 
gas -phase  carbon  species  increases  with  increasing  departures  from  equilib- 
rium. However,  even  for  the  relatively  large  values  of  heat  transfer  coeffi- 
cient shown  in  Figs.  11  through  13,  the  difference  between  the  two  assumptions 
is  insignificant.  It  is  concluded,  therefore,  that  for  nominal  ICBM  reentry 
situations,  the  use  of  the  frozen  assumption  will  yield  about  the  same  reces- 
sion as  the  equilibrium  assumption. 

It  should  be  noted  that  these  results  for  the  frozen  carbon  species  are  in 

2 

disagreement  with  the  calculations  of  Maurer  . That  is,  Maurer's  results 
indicate  that  use  of  the  frozen  assumption  decreases  the  value  of  B'  below  the 
value  obtained  for  the  equilibrium  assumption,  whereas  the  present  results 
indicate  larger  B*  at  low  temperatures  and  smaller  B'  at  high  temperatures. 
Maurer  does  not  describe  how  his  frozen  results  were  obtained,  and  this 
discrepancy  is  currently  unresolved. 

E.  CARBON  MELT  LAYER  CALCULATIONS 

The  mass  transfer  calculations  for  vaporization  from  a liquid  carbon 
surface  are  presented  herein.  As  indicated  earlier,  the  critical  assumption 
used  in  these  calculations  is  that  the  liquid  carbon  layer  is  thin  enough  so  that 
it  adheres  to  the  surface,  does  not  flow,  and  is  not  removed  by  aerodynamic 
forces.  The  calculations  were  carried  out  using  the  analysis  of  Subsection 
11-A  with  the  liquid  carbon  vapor  pressure  law  of  the  Appendix. 

Since  the  liquid  carbon  vapor  pressure  is  always  smaller  than  the  solid 
carbon  vapor  pressure,  the  mass  transfer  from  the  liquid  surface  should  be 
less  than  that  from  the  solid  surface  at  the  same  conditions.  This  is  evident 
in  P'igs,  14  and  15  where  the  equilibrium  values  of  B'  are  presented  for  two 
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Fig.  15.  Equilibrium  Vaporization  from  a 
Carbon  Melt  Layer  (P=  100  atm) 
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pressures  and  several  melt  temperatures.  The  nonequilibrium  calculations 
show  the  same  trend  as  seen  in  Fig.  16. 

Note  that  implicit  in  these  results  is  the  assumption  that  the  triple 
point  pressure  is  less  than  the  indicated  pressure  for  which  the  calculations 
were  made.  That  is,  liquid  carbon  cannot  exist  in  equilibrium  at  pressures 
lower  than  the  triple  point  pressure.  Thus,  for  example,  if  the  triple  point 
pressure  is  of  the  order  of  100  atm,  then  the  equilibrium  calculations  at  one 
atm  (Fig.  14)  are  meaningless.  The  correct  value  of  the  triple  point  pres- 
sure is  currently  in  dispute,  with  quoted  values  ranging  from  less  than  one 
19  20 

atm  (Whittaker  ) to  values  over  100  atm  (Gokcen  ).  Therefore,  calcula- 
tions are  presented  herein  for  liquid  carbon  over  this  range  of  pressures. 

Although  the  results  indicate  large  differences  in  recession  rate, 
dependent  upon  the  surface  phase  state,  the  solution  of  the  surface  energy 
balance  suppresses  the  mass  loss  difference  while  yielding  substantial  dif- 
ferences in  the  predicted  surface  temperature.  This  is  indicated  in  Figs.  17 
and  18  for  equilibrium  sublimation,  using  the  steady-state  conduction  model 
in  the  surface  energy  balance.  It  is  obvious  that  the  difference  in  mass  loss 
is  insignificant,  whereas  the  temperature  difference  is  about  1000°R.  There- 
fore, if  thermochemical  mass  loss  were  the  only  ablation  mechanism,  the 
surface  phase  question  would  be  somewhat  academic.  However,  the  erosion 
mass  loss  mechanisms  are  dependent  upon  surface  temperature  which,  as 
indicated  by  Fig.  18,  is  influenced  by  the  surface  state. 


iq 

Whitaker,  A.G.  and  P.  L.  Kintner,  "Laser  Heating  Studies  of  Carbon 
Melting  and  Vaporization,  " American  Ceramic  Society  27th  Pacific  Coast 
Regional  Meeting,  North  Hollywood,  Calif.,  October  1974. 

Gokcen,  N.A.,  et  al.,  "Determination  of  Graphite/Liquid/ Vapor  Triple 
Point  by  Laser  Heating,  " High  Temperature  Science,  8,  pp.  81-97,  1976. 


-58- 


in.  / sec 


Fig.  17.  Comparison  of  Solid  and  Liquid  Carbon 
Ablation  Models  for  Surface  Recession 
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Also  shown  in  Fig.  17  is  the  upper-bound  recession  rate  associated 
with  mechanical  removal  of  the  melt  layer  before  it  vaporizes.  The  steady- 
state  surface  energy  balance  yields  the  following  expression  for  the  mechani- 
cal removal  model; 
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(48) 


That  is,  the  surface  is  presumed  to  be  at  the  melt  temperature,  such  that  the 

mass  removed  absorbs  only  the  energy  required  to  raise  it  from  to  T 

B m 

plus  the  heat  of  fusion.  Since  there  is  no  vaporization  of  the  liquid,  there  is 

no  blowing  correction  to  the  aerodynamic  heat  flux,  and  the  gas  adjacent  to 

the  surface  is  assumed  to  be  air  with  enthalpy  evaluated  at  the  surface 

temperature  (i.e.,  T ). 

m 

Therefore,  Fig.  17  contains  both  the  upper-  and  lower-bound  recession 
rate  models  from  a liquid  carbon  surface.  It  is  apparent  that  the  upper-bound 
recession  rate  is  only  about  2.  5 times  greater  than  the  lower-bound  value. 

This  relatively  small  difference  (compared  to  other  ablative  materials)  stems 
from  the  large  heat  of  fusion  (see  Appendix)  and  the  high  melt  temperature, 
which  results  in  a large  effective  heat  of  ablation.  For  the  melt  temperature 
of  6840°R  and  heat  of  fu  sion  of  25  kcal/mol,  the  denominator  of  Eq.  (48)  yields 

H,(T  ) - = 6590  Btu/lb  (49) 

Lj  m o 

It  was  speculated  in  Section  II- D that  the  large  recession  rates  predicted 
by  Kratsch  were  due  to  mechanical  removal  of  the  melt  layer.  The  results 
obtained  from  Eq.  (48)  indicate  much  smaller  recession  rates,  indicating  that 
mechanical  removal  is  not  the  explanation.  In  any  event,  it  is  obvious  that 
the  present  results  are  in  total  disagreement  with  those  of  Kratsch. 
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The  steady-state  melt  layer  thickness  obtained  from  Eq.  (44)  is  shown 

in  Fig.  19  for  two  different  values  of  recovery  enthalpy.  The  liquid  thermal 

diffusivity  was  assumed  to  be  an  order  of  magnitude  smaller  than  the  diffu- 

2 1 

sivity  of  solid  graphite  . Since  the  melt  layer  thickness  is  linear  with  the 

liquid  diffusivity,  the  thickness  for  other  values  of  a,  may  be  obtained  by 

1-1 

scaling  the  values  of  Fig.  19. 

The  results  indicate  that  the  melt  layer  thickness  for  a flight  environ- 
ment will  range  from  one  mil  to  10  mile  depending  upon  the  severity  of  the 
environment.  Although  these  values  are  quite  small  relative  to  a typical 
nosetip  radius,  it  cannot  be  stated  that  they  are  sufficiently  small  to  resist 
flowing  of  the  liquid  carbon.  Resolution  of  this  question  would  require  an 
analysis  of  the  dynamics  of  the  liquid  layer  and  is  beyond  the  scope  of  this 
report.  The  calculation  for  a low  enthalpy  flow  (typical  of  a high-pressure 
arc  jet)  yields  melt  layer  thicknesses  less  than  half  a mil.  The  greatly 
reduced  thickness  for  the  low  enthalpy  case  is  due  to  a much  lower  surface 
temperature  (=s7000^R)  than  that  of  Fig.  18.  These  small  thicknesses  may  j 

explain  why  molten  carbon  is  not  readily  observable  in  high-pressure  arc  jet  j 

facilities.  ! 

I 

F.  STEADY-STATE  CONDUCTION  ASSUMPTION  i 

As  indicated  in  Subsection  II- E,  the  use  of  the  steady- state  conduction  j 

assumption  can  lead  to  substantial  errors  in  the  predicted  surface  mass  loss. 

In  this  subsection,  calculations  are  presented  for  both  the  steady-state  con- 
duction model  and  the  transient  expression  of  Eq.  (41).  Nominal  graphite 
properties  were  used,  and  the  thermal  diffusivity  was  taken  as  0.  0001  ft  /sec. 

The  surface  ene  rgy  balance  was  solved  for  several  values  of  heat  flux  using 
the  equilibrium  sublimation  assumption.  The  results  are  shown  in  Fig.  20, 
where  it  is  apparent  that  the  use  of  the  steady- state  conduction  model  over- 
predicts  the  recession  by  as  much  as  a factor  of  two  at  early  times.  As 

^^White,  H.  M.,  "Estimation  of  Some  Carbon  Properties  at  the  Triple  Point,  " 

Report  No.  TOR -66 9(S68 1 1 - 20) - 12,  The  Aerospace  Corporation,  San 
Bernardino,  Calif.,  June  1966. 
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indicated  by  the  criteria  of  Eq.  (42),  the  steady- state  approximation  becomes 
increasingly  valid  as  time  and  recession  rate  increase. 


G.  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  equilibrium  and  nonequilibrium  ablation  models  are  compared,  in 

this  subsection,  with  the  available  carbon  sublimation  data.  Since  the  main 

cotjcern  of  this  report  is  with  convective  heating  in  an  air  environment, 

attention  is  restricted  to  data  which  fall  into  this  category.  The  best  such 

data  has  been  obtained  at  NASA-AMES  in  a low-pressure  arc  jet  facility, 

22  23 

first  by  Lundell  and  later  by  Maurer  . The  main  difference  between  the 
two  sets  of  data  is  that  the  convective  heat  flux  for  Maurer's  tests  was  aug- 
mented by  a radiative  heat  flux  to  obtain  higher  surface  temperatures.  The 
only  correction  made  to  the  data  was  to  add  100°K  to  the  tabulated  surface 
brightness  temperature  of  Maurer  to  obtain  the  true  temperatures.  The 
blowing  law  of  Eq.  (13)  was  used  to  cast  the  data  in  terms  of  B'. 

The  equilibrium  and  nonequilibrium  calculations  are  compared  with  the 
data  in  Figs.  21  through  23.  The  nonequilibrium  calculations  use  the  slow 
nitridization  assumption  and  the  Dolton  vaporization  coefficients.  The  Lundell 
data  (Fig.  21)  is  in  good  agreement  with  the  equilibrium  theory.  This  close 
agreement  obviously  precludes  the  necessity  of  invoking  mechanical  erosion 
arguments  to  match  the  data,  (This  does  not,  of  course,  explain  the  presence 
of  the  pjarticulate  matter  observed  by  Lundell.) 

In  contrast  to  the  Lundell  data,  the  Maurer  data  is  seen  to  be  in  closer 
agreement  with  the  nonequilibrium  calculations.  The  low-pressure  data  of 
Fig.  23  lies  below  the  lower  bound  melt  temperature  of  6840°R  and  is  in 
excellent  agreement  with  the  nonequilibrium  calculation  for  solid  carbon. 


Lundell,  J.H.  and  R.  R.  Dickey 
Temperatures,  " AIAA  Journal, 


"Ablation  of  ATJ  Graphite  at  High 
11(2),  pp.  216-222,  January  1973. 


Maurer,  R.  E.,  et  al.,  "Graphite  Sublimation  Under  Low  and  High  Con 


vective  Mass  Transfer  Environments,  " ASME  Paper  76-ENAS-68, 
July  1976. 
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However,  the  data  of  Fig.  22  lies  above  this  value  for  melt  temperature  and 
is  in  good  agreement  with  the  nonequilibrium  calculation  for  liquid  carbon. 
While  this  may  tend  to  support  the  presence  of  liquid  carbon  (if  the  triple 
point  pressure  is  less  than  one  atm),  this  conclusion  can  hardly  be  regarded 
as  other  than  very  tentative.  No  evidence  of  liquid  carbon  was  seen  during 
the  tests  (the  predicted  steady-state  melt  layer  thickness  was  0.003  in.  ) and, 
in  particular,  the  quoted  accuracy  of  the  surface  temperature  measurements 
(±2  percent)  is  sufficiently  large  to  move  the  data  to  the  solid  nonequilibrium 
calculation  of  Fig.  22. 

Even  allowing  for  the  uncertainty  in  surface  temperature,  the  Maurer 
data  shows  definite  nonequilibrium  effects;  it  is  concluded  that  the  nonequilib- 
rium calculation  using  the  Dolton  vaporization  coefficients  best  represents 
the  data. 


-70- 


IV.  SUMMARY  AND  CONCLUSIONS 


An  examination  of  the  nonequilibrium  carbon  ablation  process  has  been 
presented.  The  study  has  focused  on  the  shortcomings  of  the  available 
analyses  and  has  undertaken  a systematic  examination  of  the  assumptions 
that  are  usually  invoked  in  the  nonequilibrium  model  for  carbon  ablation. 

The  major  areas  of  investigation  included  the  influence  of  the  assumed 
chemical  state  (i.e.,  frozen  or  equilibrium  reactions),  ablation  from  a mol- 
ten carbon  surface,  and  the  computational  considerations  associated  with 
solving  the  system  of  equations. 

The  assumptions  of  fast  and  slow  heterogeneous  nitridization  reactions 
were  examined  in  conjunction  with  an  equilibrium  gas  phase.  The  results 
indicated  significant  recession  rate  differences  only  for  extreme  departures 
from  equilibrium.  It  was  concluded  that,  for  the  nominal  range  of  pressure 
and  heat  flux  for  ground  test  and  flight  environments,  the  differences  are 
small  (on  the  order  of  10  percent).  The  assumption  that  the  carbon  sublima- 
tion species  are  frozen  in  the  gas  phase  was  also  examined  and  found  to  exert 
no  significant  influence  on  the  recession  rate. 

Although  the  various  reaction  rate  assumptions  were  found  to  exert 
only  small  influences  on  the  recession  rate,  the  only  self-consistent  model 
was  found  to  be  the  slow  nitridization  assumption  with  an  equilibrium  gas 
phase.  Comparisons  of  the  slow  nitridization  model  with  experimental  data 
show  reasonable  agreement  (using  Dolton  vaporization  coefficients),  and  this 
model  is  recommended  for  adoption. 

The  sublimation  from  a liquid  carbon  surface  was  handled  by  assuming 
that  the  melt  layer  was  sufficiently  thin  to  adhere  to  the  surface  without 
mechanical  removal.  The  liquid  carbon  assumption  was  found  to  result  in 
a substantial  reduction  of  the  predicted  B',  as  compared  to  the  value  ob- 
tained for  vaporization  from  a solid  surface.  Application  of  the  carbon  melt 
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layer  calculation  to  flight  environments  yielded  slightly  lower  recession 
rates  than  the  solid  carbon  assumption  but  with  an  increase  in  surface 
temperature  of  about  1000°R.  It  is  hoped  that  these  results  will  provide 
some  impetus  for  the  experimental  investigation  of  the  liquid  layer  hypothe- 
sis in  ground  test  arc  jet  facilities. 

A unified  analysis  of  the  nonequilibrium  carbon  ablation  process  has 
been  presented  by  treating  the  oxidation  and  sublimation  regimes  simul- 
taneously. A computationally  efficient  solution  procedure,  requiring  only  a 
single  iteration  loop,  is  utilized  to  solve  the  system  of  equations.  Conver- 
gence of  the  iteration  scheme  is  sufficiently  fast  to  utilize  the  method  as  a 
subroutine  package  for  inclusion  into  thermodynamic  material  response  codes. 

Finally,  some  comments  concerning  recommendations  for  future  abla- 
tion tests  are  in  order.  The  analysis  of  this  report  has  treated  the  vapori- 
zation coefficients  as  known  quantities.  The  results  of  this  investigation 
indicate  that  the  uncertainties  in  the  vaporization  coefficients  exert  sub- 
stantially more  influence  on  the  ablation  predictions  than  do  the  uncertainties 
in  the  gas/solid  reaction  rates.  In  order  to  experimentally  investigate  the 
vaporization  process  in  the  absence  of  chemical  reactions,  it  is  recommended 
that  a series  of  arc  jet  ablation  tests  be  carried  out  in  an  inert  atmosphere. 
Both  ATJ-S  and  the  current  carbon-carbon  composite  materials  should  be 
tested.  In  addition,  since  no  carbon-carbon  ablation  data  exists  in  the  sub- 
limation regime,  arc  jet  tests  are  also  required  for  this  material  in  air. 

Further  areas  of  investigation  should  include  experiments  to  determine 
the  C^N  heat  of  formation  and  reaction  rates  for  heterogeneous  carbon- 
nitrogen  reactions,  for  both  graphite  and  carbon-carbon. 


i 

I -72- 


i 


REFERENCES 


1.  Crowell,  P.G.,  "Analytic  Solutions  for  Carbon  Sublimation  in  Atmos- 
pheres of  Elemental  Carbon,  Hydrogen,  Nitrogen  and  Oxygen,  " 

Report  No.  TR-0078(3550-15)-2,  The  Aerospace  Corporation,  El 
Segundo,  Calif.,  December  1977. 

2.  Maurer,  R.  E.,  "An  Evaluation  of  Some  Uncertainties  in  Carbon  Sub- 
limation Kinetics,"  Aerotherm  TM-73-27,  January  1973. 

3.  Kratsch,  K.  M.,  "Graphite  Fusion  and  Vaporization,  " ASME/AIAA  10th 
Structures,  Structural  Dynamics  and  Materials  Conference,  April  1969. 

4.  Dolton,  T.A.,  R.  E.  Maurer,  and  H.  E.  Goldstein,  "Thermodynamic 
Performance  of  Carbon  in  Hyperthermal  Environments,  " AIAA  3rd 
Thermophysics  Conference,  Paper  No,  68-754,  June  1968. 

5.  Baker,  R.  L..,  "Graphite  Sublimation  Chemistry  Nonequilibrium  Effects,  " 
AIAA  Journal.  ^(10),  October  1977. 

6.  Ziering,  M.  B.  and  V.  DiCristina,  "Thermomechanical  Erosion  of 
Ablative  Plastic  Composites,  * AIAA  7th  Thermophysics  Conference 
April  1972,  Paper  No.  72-299. 

7.  JANNAF  Thermochemical  Tables,  National  Bureau  of  Standards 
NBS-37.  June  1971. 

8.  Kennard,  E.H.,  "Kinetic  Theory  of  Gases,  " McGraw-Hill  Book  Co. 

Inc.,  New  York,  1938. 

9.  Goldstein,  H.  W.,  "The  Reaction  of  Active  Nitrogen  with  Graphite,  " 
Journal  of  Physical  Chemistry.  ^(1),  January  1964. 

10.  Dorrance,  W.H.,  "Viscous  Hypersonic  Flow,  " McGraw-Hill  Book  Co. 

Inc.,  New  York,  1962.  ' 

11.  Carslaw,  H.S.  and  J.C.  Jaeger,  Conduction  of  Heat  in  Solids.  Oxford 

University  Press,  1959,  p.  388.  ^ 

12.  Hansen.  C.F.  and  W.  E.  Pearson,  "A  Quantum  Model  for  Bending 

Vibrations  and  Thermodynamic  Properties  of  C,,  " Canadian  Journal  of 
Physics,  51,  1973.  ^ 

13.  Pearson,  W.E,  and  W.  C.  Davy,  "New  Thermodynamic  Functions  for 
the  C3  Molecule,  " AIAA  Journal.  1_1(8),  August  1973. 


-73- 


REFERENCES  (Continued) 


14.  Strauss,  H.  L.  and  E.  Thiele,  "Thermodynamics  of  C3  General  Methods 
for  Nonrigid  Molecules  at  High  Temperature,  " Journal  of  Chemical 
Physics,  ^(7),  April  1967. 

15.  Wachi,  F.M.  and  D.  E.  Gilmartin,  "Heat  of  Formation  and  Entropy  of 
C3  Molecule,  " Report  No.  TR-0172(2250-40)-7,  The  Aerospace  Corpo- 
ration, El  Segundo,  Calif.,  May  1972. 

16.  Meyer,  J.A.  and  D.  W.  Setser,  "Chemiluminescence  from  the  Reaction 

of  Oxygen  Atoms  with  Dicyanoacetylene,  " Journal  of  Physical  Chemistry. 
74,  1970,  p.  3452.  

17.  Hirth,  J.  P.,  "Evaporation  and  Sublimation  Mechanisms,  " The  Charac- 
terization of  High  Temperature  Vapors,  J.  L.  Margrave  (ed.),  John 
Wiley  & Sons,  Inc.,  New  York,  1967, 

18.  Bernstein,  H.  and  J.  R.  Baron,  "Stability  Characterization  of  Refractory 
Materials  Under  High  Velocity  Atmospheric  Flight  Conditions,  Part  IV, 
Volume  II,  Calculation  of  the  General  Surface  Reaction  Problem,  " 

AFM L TR -69- 84,  December  1969. 

19.  Whitaker,  A.G.  and  P.  L.  Kintner,  "Laser  Heating  Studies  of  Carbon 
Melting  and  Vaporization,  " American  Ceramic  Society  27th  Pacific 
Coast  Regional  Meeting,  North  Hollywood,  Calif.,  October  1974. 

20.  Gokcen,  N.A.,  et  al.,  "Determination  of  Graphite/ Liquid/ Vapor  Triple 
Point  by  Laser  Heating,  " High  Temperature  Science,  8,  pp.  81-97. 

1976. 

21.  White,  H.M,,  "Estimation  of  Some  Carbon  Properties  at  the  Triple 
Point,  " Report  No.  TOR -669(S681 1 -20)- 12,  The  Aerospace  Corporation, 
San  Bernardino,  Calif.,  June  1966. 

22.  Lundell,  J.H.  and  R.  R.  Dickey,  "Ablation  of  ATJ  Graphite  at  High 
Temperatures,  " AIAA  Journal.  H^(2),  pp.  216-222,  January  1973. 

23.  Maurer,  R.  E.,  et  al.,  "Graphite  Sublimation  Under  Low  and  High 
Convective  Mass  Transfer  Environments,  " ASME  Paper  76-ENAS-68 
July  1976. 

24.  Powars,  C.A.  and  R,  M.  Kendall,  "User's  Manual- -Aerotherm  Chemi- 
cal Equilibrium  (ACE)  Computer  Program,  " Aerotherm  Corporation. 

May  1969. 


-74- 


APPENDIX  A 
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VAPOR  PRESSURE  OF  CARBON  SUBLIMATION  SPECIES 


As  discussed  in  the  main  text  (Section  II),  the  vapor  pressure  of  each 
gaseous  carbon  species  is  assumed  to  be  a known  function  of  temperature. 
These  vapor  pressures  are  presented  in  this  Appendix  for  vaporization  from 
both  solid  and  liquid  carbon  surfaces.  The  solid  vaporization  case  is  con- 
sidered first. 

A.  VAPORIZATION  FROM  SOLID  CARBON 

The  gaseous  carbon  species  in  equilibrium  with  the  solid  carbon  surface 
are  denoted  by  C^.  The  vaporization  reaction  is  given  as 


'solid 


(A-1) 


For  this  reaction,  the  third  law  gives 


in  P, 


- E„(Ci)  /f“  - 


r»-H» 


(A-2) 


solid 


where  Hq{C^)  is  the  heat  of  formation  of  Cj^  at  0°K,  and  the  heat  of  formation  of 
solid  carbon  has  been  set  equal  to  zero  by  convention. 

The  free -energy  function  and  partition  function  for  the  gaseous  species 
are  related  as  follows: 


The  free-energy  function  is  related  to  the  specific  heat  by 


f 

i 

i 

V 


^ ■ ^298 


• T / S"  - / 


^ - =298 


(A -4) 


298 


298 


Note  that  the  JANNAF  thermochemical  data  uses  a base  temperature  of 
2 98°K.  The  free-energy  functions  at  0°K  and  298°K  are  related  by 


r»-H» 


F° 

^ ■ ^298 


"O  ■ 298 


(A-5) 


The  free-energy  function  referenced  to  298°K  and  the  difference  in 
enthalpies  between  0°K  and  298°K  are  both  tabulated  in  JANNAF. 

For  solid  carbon,  the  JANNAF  data  yields 


- hJ 


F° 

■ 298  252 

rm  ' rrs 


(A-6) 


The  curve  fits  for  specific  heat  in  the  ACE  code  were  used  to  evaluate 
the  free-energy  function  for  solid  carbon  by  using  Eq.  (A-4).  The  final 
result  is 


F"  - H 


298 


= A(1  - inT)  - Y BT  - 


2T 


(A -7) 


T < 3000°K: 

A = 5.  861 
B = 0.  954  X lO’"^ 
C = -0. 7666  X 10^ 
D = 4323. 3 
E = 35.  235 
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T ^ 3000°K: 


A *4.85 
B = 0. 2916  X 10’^ 

C = 0. 3072  X 10^ 

D = 893.7 
E = 27.  516 

where  all  temperatures  are  in  °K,  and  the  units  of  the  free-energy  functions 
are  cal/mol-°K. 

The  use  of  Eqs.  (A-6)  and  (A-7)  in  Eq.  (A-2)  provides  all  the  necessary 
information  about  the  solid  phase. 

The  partition  functions  were  used  to  obtain  the  free-energy  functions  for 
the  gaseous  carbon  species  (with  the  exception  of  C3).  The  free  energy 
function  was  obtained  from  the  curve  fit  by  Pearson^^. 

The  enthalpy  of  solid  carbon  is  required  in  the  surface  energy  balance 
and  is  given  by 

H°(S)  ^ AT  BT^  - £ + 252  - D cal/mol  (A-8) 

where  (S)  denotes  solid  properties. 

B.  VAPORIZATION  FROM  LIQUID  CARBON 


For  vaporization  from  a liquid  carbon  surface,  the  third  law  yields 


The  liquid  carbon  free  energy  function  is  given  by 


liquid 


k f 


C„(L) 


J Cp(L)dT  - J 

T T 

m m 


dT  - (L) 


(A-10) 


The  entropy  and  enthalpy  of  the  liquid,  at  the  melt  temperature,  are  related  to 
the  corresponding  solid  values  through  the  heat  of  fusion. 


(L)  = (S)  + 

m m 


S°  (L)  = S®  (S)  + 

m m m 


(A-11) 


The  enthalpy  and  entropy  of  the  solid,  at  the  melt  temperature,  are 
obtained  from  the  JANNAF  data. 


The  additional  data  required  to  evaluate  these  expressions  is  the  melt 


temperature,  the  specific  heat  of  liquid  carbon,  and  the  heat  of  fusion  at  the 


melt  temperature.  The  following  values  were  obtained  from  White  : 


= 25  Kcal/mol 
F 


Cp(L)  = 7.44  cal/mol-  K 


(A-12) 


For  a constant  value  of  the  liquid  carbon  specific  heat,  the  free-energy 
function  of  Eq.  (A-10)  becomes 


fO  -h£ 


«T  \ 

m I 

/liquid 


T 1 n 

= C (L)  l-^+/n^  - S!J,  (L)  (A-13) 

^ L J m 
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The  lower-bound  melt  temperature  of  3800°K,  obtained  by  Whittaker 
was  used  in  most  of  the  calculations.  The  JANNAF  values  for  solid  carbon,  at 
this  temperature,  are  given  as 


T = 3800“K 
m 


(S)  = 19559  cal/mol 


(S)  = 13.  575  cal/mol-°K 


(A-14) 


The  enthalpy  of  the  liquid  carbon  surface  is  required  to  solve  the  sur- 
face energy  balance  and  is  given  as 


H°(L) 


= (S)  + + /*  C (L)  dT 


(A-15) 


All  the  calculations  of  this  report  used  a constant  value  of  C (L),  for  which 
the  integral  of  Eq.  (A-15)  is  easily  evaluated. 

Finally,  it  is  instructive  to  examine  the  ratio  of  liquid  to  solid  vapor 
pressures  for  temperatures  greater  than  the  melt  temperature.  The  third 
law  relationships  of  Eqs.  (A-1)  and  (A-9)  are  combined  to  yield 


r 0 0 

^ - i F^(S) 

(S)  RT  ■ RT 

1 


(A-16) 


Substitution  of  the  free-energy  functions  into  this  expression  yields 

5 j + T J - '=p(S'l  ‘‘T 


C„(L)  - C (S) 


(A-17) 


r 


It  is  apparent  that  the  second  integral  is  larger  in  magnitude  than  the 

first  integral  provided  that  the  liquid  specific  heat  exceeds  the  solid  specific 

heat.  This  is  clearly  the  case  for  the  value  of  Eq.  (A- 12)  (as  compared  with 

the  tabulated  JANNAF  values).  Thus,  for  T >T  it  is  obvious  that  the  right- 

m ® 

hand  side  of  Eq.  (A- 17)  is  negative.  Therefore,  the  vapor  pressure  of  carbon 
I species  C^,  in  equilibrium  with  a liquid  carbon  surface,  is  less  than  the  vapor 

I pressure  of  the  same  species  in  equilibrium  with  a solid  carbon  surface  at  a 

temperature  greater  than  the  melt  temperature.  That  is 

P (L)<P  (S)  (A- 18) 

i ^i 

I 
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APPENDIX  B 


SELECTION  OF  CHEMICAL  SPECIES 

At  the  suggestion  of  the  review  committee,  additional  discussion  and 
justification  for  the  choice  of  included  species  is  provided  in  this  appendix. 

The  choice  of  allowable  species  was  constrained  by  accuracy  and  com- 
putational considerations.  Thus,  all  of  the  dominant  species  identified  by 
previous  investigators  were  automatically  included.  Only  two  minor  species 
were  included  (atomic  oxygen  and  nitrogen)  because  they  are  potentially 
important  with  respect  to  the  heterogeneous  oxidation  and  nitridization  reac- 
tions, and  their  inclusion  does  not  impact  the  solution  procedure,  ! 

The  inclusion  of  nitrogen-oxygen  species  into  the  formulation  intro- 
duces additional  iteration  loops  into  the  solution  procedure,  is  computationally 
expensive,  and  is  not  warranted  in  terms  of  accuracy.  At  high  temperatures, 
the  nitrogen  must  compete  with  the  carbon  for  the  available  oxygen,  and  since 
oxygen  reacts  far  more  readily  with  carbon,  there  is  very  little  excess  oxygen 
available  for  the  formation  of  nitrogen -oxygen  compounds.  At  low  tempera- 
tures, where  excess  oxygen  is  available,  the  relevant  equilibrium  constants 
are  so  small  that  nitrogen-oxygen  compounds  exist  only  as  trace  species. 

These  intuitive  considerations  were  checked  by  running  some  exact  equilibrium 
24 

cases  with  the  ACE  code  over  the  range  of  pressures  and  temperatures  of 
interest.  The  results  verify  that  the  nitrogen-oxygen  compounds  ai^  insigni- 
ficant. Thus,  none  of  these  species  are  included  in  the  formulation. 

The  ACE  results  also  indicur-;  that  the  only  significant  carbon -oxygen 
species  are  CO  and  CO^  and  these  are  the  only  such  species  included. 


24 

Powars,  C.A.  and  R.M.  Kendall,  "User's  Manual- -Ae rotherm  Chemical 
Equilibrium  (ACE)  Computer  Program,  " Ae  rotherm  Corporation, 

May  1969. 


-81- 


The  only  carbon  sublimation  species  included  were  those  for  which 
reliable  spectroscopic  data  exists  (i.e.,  C^-C^).  Although  some  investi- 
gations have  included  carbon  species  up  to  their  properties  were 

obtained  from  theoretical  models  utilizing  assumption  of  unknown  accuracy. 
Although  the  inclusion  of  the  higher  order  carbon  species  does  not  impose 
any  additional  computational  difficulty,  it  was  felt  that  the  uncertainty  asso- 
ciated with  their  thermochemical  properties  data  justified  their  deletion. 


NOMENCLATURE 


I 


f 


[ 


B' 


H 


R 


H, 


H 


K, 


K. 

1 

K. 


M 

M, 


M. 

1 


M 


C. 


m 


m. 

1 


P 

P. 


Pc. 

1 


R 

S 

T 

t 


mass  transfer  parameter  [defined  in  Eq.  (3)] 

enthalpy  of  solid  carbon  evaluated  at  some  specified  bulk 
temperature 

enthalpy  of  liquid  carbon  at  temperature,  T 
recovery  enthalpy 

enthalpy  of  carbon  surface  phase  (either  liquid  or  solid) 

enthalpy  of  gas  mixture  adjacent  to  surface 

pressure  equilibrium  constant 

mass  fraction  of  species  i 

mass  fraction  of  element  i 

I 

molecular  weight  of  gas  mixture 
molecular  weight  of  inert  species 
molecular  weight  of  species  i 

molecular  weight  of  carbon  sublimation  species  C^^ 
total  carbon  mass  loss 

mass  loss  due  to  heterogeneous  reaction  of  carbon  with  species  i 

surface  pressure 

partial  pressure  of  species  i 

equilibrium  vapor  pressure  of  carbon  sublimation  species 

cold  wall,  nonblowing  heat  flux 

universal  gas  constant 

surface  recession  rate 

temperature 

time 
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T 


NOMENCLATURE  (Continued) 


6 

P U c„ 

e H 

P u C„ 
'^e  e H 


0 


a 


vaporization  coefficients  for  carbon  sublimation  8p>ecie8 
carbon  melt  layer  thickness 

cold  wall,  blowing  heat  flux  parameter  = q/H^ 
cold  wall,  nonblowing  heat  flux  parameter  = q^/Hp 
thermal  diffusivity 


Subscripts 

e 

L 

w 


boundary  layer  edge  values 

liquid  carbon 

wall  (surface)  values 


Superscripts 

* denotes  a heterogeneous  carbon-nitrogen  reaction 

~ denotes  an  element 


I 
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